ON THE DISTRIBUTION OF IMAGINARY PARTS OF ZEROS OF THE 

RIEMANN ZETA FUNCTION 



KEVIN FORD AND ALEXANDRU ZAHARESCU 

Abstract. We investigate the distribution of the fractional parts of cry, where a is a fixed 
non-zero real number and 7 runs over the imaginary parts of the non-trivial zeros of the 
Riemann zeta function. 



1. Introduction 

There is an intimate connection between the distribution of the nontrivial zeros of the 
Riemann zeta function ((s) and the distribution of prime numbers. Critical to many prime 
number problems is the horizontal distribution of zeros; here the Riemann Hypothesis (RH) 
asserts that the zeros all have real part |. There is also much interest in studying the 
distribution of the imaginary parts of the zeros (the vertical distribution). For example, one 
expects that their consecutive spacings follow the GUE distribution from random matrix 
theory. Originally discovered by Montgomery [12J, who studied the pair correlation of zeros 
of the zeta function, this phenomenon has been investigated, for higher correlations and 
also for more general L— functions, by a number of authors, including Odlyzko ^Hj, Hejhal 
jO], Rudnick and Sarnak J7|, Katz and SarnakJSJ, Murty and Perelli [13], and Murty and 
Zaharescu |14j . 

Let {y} denote the fractional part of y, which can be interpreted as the image of y in the 
torus T = R/Z. In this paper we look at the distribution of {cry} where a is a fixed nonzero 
real number and 7 runs over the imaginary parts of the zeros of ((s). The starting point is 
an old formula of Landau ^T], which states that for each fixed x > 1, 

(1.1) x p = -^A(x) + 0(\o g T). 

0<7<T 

Here the sum is over zeros p = (3 + ry of (, A (2) is the von Mango ldt function for integral 
x > 1 and A (a;) = for non- integral x > 1. Put x = e 27rja into where a > and j is a 

positive integer. If RH holds, then (jl.lj) implies that 

(1.2) e 2mja7 = 0(T), 

0<7<T 
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while there are ^> TlogT summands on the left side (see (jl.6j) below). Thus, as Rademacher 
observed in 1956 ^B], RH plus Weyl's criterion (PU], Ch. 1, Theorem 2.1) implies that the 
sequence {07} is uniformly distributed in [0, 1). In 1975, Hlawka 7J proved this conclusion 
unconditionally. His result depended on a version of ()1.1|) which is uniform in x and results 
about the density of zeros lying off the line (3 = 5 ( see Lemmas 1,2 below). Hlawka also 
showed via the Erdos-Turan inequality (jlUj. Ch. 2, Theorem 2.5) that the discrepancy 
of the set {{07} : < 7 < T} is O(j^y) assuming RH. More recently Fujii 3] showed 
unconditionally that this discrepancy is 0( 1 ° 1 ^ 1 °|, T ). 

Ramemacher also noted in ^1 that if a = ° sp , where p is a prime and k e N (i.e. 
x = p k so A(x) = logp), by ()1.1|) there should be a "predominance of terms which fulfill 
({cry} — 1/2 1 < 1/4". We will give a very precise meaning to this statement, and also show 
that Hlawka's discrepancy result is best possible in a certain sense. We will see that there is 
a certain limiting measure fi a that one can naturally associate to each a > 0, which explains 
among other things the behavior of the above discrepancy as T — > 00. 

For any real numbers a, T > consider the measure \x a% T defined on T by 

. . 1 r N ( T ) 

(!-3) fJ, a>T = — £{cry} j—V- 

0<7<T 

Here \i denotes the Haar measure on T, <5{ a7 } is a unit point delta mass at {07} and N(T) 
is the number of zeros p — (3 + ry of (, counted with multiplicity, with < j3 < 1 and 
< 7 < T. To the measure /z Qj t we associate the function 



;i.4) M(y;T) = M a (y;T):=^ T ([0,y)) = ± ^ 1 



T 

0<7<T 
{a-y}<y 



N(T) 



In particular, M(0;T) = M(l; T) = 0. Let 



(1-5) D* a (T) = D*(T) = sup \M{y-T)\ 

0<y<l ^ {-L ) 

denote the discrepancy of the set {{07} : < 7 < T}. 

For fixed a, the above measures fM a) T, seen as continuous linear functionals on C(T), form 
an unbounded set. In fact, the total variation of \x a ^ grows like logT as T — > 00, since (|18j. 
Theorem 9.4) 

(1.6) iV(T) = ^log^-^ + 0(logT) (T>1). 

We will see, however, that for a large class of functions h : T — ► C (including all absolutely 
continuous functions on RH), lirnr-Kx, L hd^x exists. Moreover, there is a unique Borel 
measure \i a on T which is absolutely continuous with respect to the Haar measure, such that 
for such h one has 

(1.7) lim / hd^ at T = / hd\i a . 

T— >oo /• / 
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Let g a : [0, 1) — > M denote, via the above identification between T and [0, 1), the density of 
this measure. Note that fi a> T(T) = for any T, hence 



;i-8) / g a (t)dt = 

o 



We have identified the density function g a , and we found that 

1.9) g a (t) = (tet) 



provided a is not a rational multiple of a number of the form with p prime. If a 



a logp 
2irq 

for some prime number p and positive integers a, q with (a, q) = 1, then for any t G [0, 1), 



9a (t) = -^l^Y^/^ntyk 
(p% cos2nqt — 1) logp 



(1.10) " 



As a function of t, g a (t) attains its global minimum at each of the points t = |, k 



7r(p a — 2p2 cos 2irqt + 1) 

; i1 

0, 1, • • ■ , q — 1, the minimum being 

/k\ logp 

1.H) 9a(-)=- , , u <0- 

\qy 7r(]?2 — 1) 

In particular, this shows that there is a shortage of zeros of ((s) with imaginary parts 7 such 
that {cry} = {3^7} is close to one of the points -, k = 0, 1, . . . , q — 1. When g = 1, this 
corresponds to Rademacher's statement mentioned above. 

We conclude this introduction by showing some histograms of M(y; T) for T = 600, 000 
(N(T) = 999508) and a few values of a. The list of zeros to this height were kindly supplied 
by Andrew Odlyzko. We partition [0, 1) into 500 subintervals of length In Figure 1 we 
plot for each subinterval I = [y, y + the value of 500(M(y + — M(y)) and also the 
graph of g a (y). 

2. Statement of results and conjectures 

To state our main theorem, we consider for any / G Lx(T) the modulus of continuity (pQ, 
Definition 1.5.1) 



u(f;5) = sup / \f(t + u)-f(t)\dt. 

\u\<5 JT 

For any real numbers a > and T > define the measure ii a ^ as in (jl.3|) , and let g a : T 
be defined by (11. 9J) or (jl.lOj) . as appropriate. 

Theorem 1. Let a > and let h : T — ► C be an absolutely continuous function with 
(2.1) uj(h';6)=o(l/\og(l/5)) (5-0+). 




Figure 1. 500(M(y + ^) - M(y)) vs. g a (y) for T = 600000. 
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Then 
(2.2) 



lim / hdjji^T = / hg a dfi. 

T — >oo 



Corollary 2. The equality Q holds for all h G C 2 (T). 
Proof. Since h G C 2 (T), 



uj(ti; 5) = sup I [ h"{y) dy dt < 5 [ \h"\ = 0{5) 



Corollary 3. Suppose a = tt 2 ° sp , where p is prime and (a, q) = 1. Then 



□ 



where 



\M{y;T)\dt > B(a) + o(l) 
J T gl dn 



B a) 



max \g' a 



Therefore, D* a (T) > B(a \ /2 ^ o{1) fo] these n. 



IobT 



Proof. Apply Theorem^ with h = g a (this may or may not be optimal). By integration by 
parts, 



9l + o(l) 



< 



M(y;T)g' a (y)dy 
\M{y;T)\dy] m^\g' a (y)\. 



□ 



Theorem 4. Let a > and suppose 

^a( T ) < i3ir- Then (D holds for all absolutely 

continuous functions h : T — > C. 

Corollary 5. On RH, ()2.2|) holds for all absolutely continuous functions /i on T. 

Proof. By Hlawka's Theorem, if RH is true then D*(T) < 1/logT. □ 

The main open problem in this line of investigation is to determine the largest "natural" 
class of functions for which ()2.2|) holds. We conjecture that (12.2}) holds for the characteristic 
function of the interval [0,y). 

Conjecture A. If a = a ^ q p , where p is prime and (a, q) = 1, then uniformly in y, 



lim M(y;T) 

1 — >oo 



g a (t) dt 



logp 

~27r 2 g 



arg (l _ p -a/2 e -2niqyy 



When a is not of this form, limr^oo M(y, T) = uniformly in y. 
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Corollary 6. Assume Conjecture A. Suppose a = a ^ q p , where p is prime and (a,q) = 1. 
As T — > oo, 

g ; ( r)^(i + o(i))^ arcs 1 in(? T :° /2) . 

7rg log i 
When a is not of this form, D* a {T) = o(l/logT). 



In the opposite direction, (|2.2|) cannot hold for all functions h which are continuous and 
differentiable on T. This is a consequence of a property of general sequences (uniformly 
distributed or not), which we state below. 

Theorem 7. Let ai, a 2 , . . . be an arbitrary sequence of numbers in T, let t be a point in T 
and let f(x) be a function decreasing monotonically to arbitrarily slowly. Then there is 
a function h, continuous and differentiable on T, and which is C°°(T\{t}), so that for an 
infinite set of n e N, 



(2.3) 



-J2h( aj )- h 

n . =1 J j 



> f(n). 



In connection with his work on /c-functions, Kaczorowski 8j was also led to study similar 
questions on the distribution of the imaginary parts of zeros of the zeta function and Dirichlet 
//-functions. In our terminology, his results concern the measures 

7>0 ' 7>0 

which essentially capture the distribution of {07} for I7 — n| y/n. Kaczorowski's Conjec- 
ture A corresponds to our Corollary |H1 and his Conjectures B, Bi and B 2 to our Conjecture 
A. 

We also mention that Fujii (|3], Theorem 5) has proved ()2.2j) in the special case h(u) = 
-B 2 (u) = u 2 — u + I (which is covered by Theorem 1) He has also proven, for a wide class of 
smooth functions /, that the sequence {/(t)} is uniformly distributed 0. 

Lastly, we mention that analogs of our results should hold for a wide variety of L-functions, 
including Dirichlet L-functions. In deriving such results, the primary tools would be a 
generalization of Lemma 1 below and zero density estimates analogous to Lemma 2 below 
which are nontrivial for certain sequences of a, T with T —>■ 00 and a — | — >• + . This will 
be addressed in a future paper. 



3. Proof of Theorem [T] 
Our principal tools are the following two lemmas. 
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Lemma 1. Let x, T > 1, and denote by n x the nearest prime power to x. Then 

, s „ Mn x ) e iTlog( > x/nx) - 1 / , 2/ mx log2T\ 

^3_]_^ \ x p = |- O I a; log (2xT) H - ) 

J-^ 2-k i\og{x/n x ) V * logx / ' 



0<7<T 

where if a; = n r the first term is — T 



A(rtg) 



2vr ' 

This is a uniform version of a theorem of Landau jTT], and the proof is nearly identical 
to the proof of Theorem 1 of Gonek ([S], §3). The only difference is in the treatment of the 
term n = n x occuring in Gonek's integral I\. That term is 

V = -^-(x/n x ) c [ T (x/n x fdt, 1 ' 



2tx J 1 log(3x) 

Gonek estimates this crudely as V <C logxmin(T, x/\x — n x \), but we are more precise. If 
\x — n x \ > 1, then V <C xlog(2x) and otherwise 

\( n ) / p iT\og(x/n x ) _ i 

V = K -^(l + 0(\x - n x \/x)) — , , . + 0(1) 

2vr 1 11 xW " V i\og{x/n x ) 1 ; 

A(n 8 ) e' rl °^) - 1 

= 0(l0g(2x)) ; r . 

1 61 " 2tt i\og{x/n x ) 

Let N(a,T) be the number of zeros p = (3 + ry of £(s) with |7| < T and /3 > a. The 
following density bound is due to Selberg (Theorem 9.19C of [TH]). 

Lemma 2. Uniformly in \ < a < 1, we have N(<r,T) < T 1 -*^- 1 /*) i og y 

We now proceed with the proof of Theorems ^ and EJ Fix a > and an absolutely 
continuous function /i : T — > C. In what follows, constants implied by the O— and <C 
—symbols may depend on a and h. Suppose either D* a {T) <C 1/logT or h satisfies (12. lj) . 
Looking at the definition of \i a ^ we see that what we need to show in order to complete the 
proof of the theorems is that 

•1 r i 



h(u)g a (u) du + o(l) 



(3.2) 1 £ M«7) ~^P- f Ku)du = I 

0<7<T ^° ^° 

as T — > 00. By treating separately the real part and the imaginary part of /i, we may assume 
in what follows that h is real. We start by approximating h by a trigonometric polynomial 
based on its Fourier series 

(3.3) h{u) = 



2nmui 



where the Fourier coefficients c m are given by 



(3.4) c m = / h(u)e- 27Tmm du. 
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Since h is absolutely continuous, we have h! G Li(T) and consequently by the Riemann- 
Lebesgue Lemma, 

(3.5) c m = o(l/m) (m — > oo). 

Let 

Kr(u)- 3 ( Sin7rjM V = 3 (Vc (1 - J+al "V 

Jl j J(2J 2 + 1) V sin™ J "J(2J 2 + 1)V^ ; 



be the Jackson kernel (pQ, Problem 1.3.9) and define 



Hj(y)= I Kj(u)h(u-y)du= ^ A^'\ ¥ ^ 

\j\<2J 



where 



A[ J) = 1, \A\ J) \ < 1, lim A { p = 1 for each fixed j. 

- 1 J— >oo J 



This kernel is chosen because it provides a very fast rate of L\ (T)-convergence of Hj to h 
(see [i] . §1.6 and Ch. 2 for more on this subject). The same is true for the convergence of 
Hj to h! . Specifically, we have (pQ, Lemma 1.5.4, Theorem 1.5.8 and Corollary 1.5.9) 

(3.6) I \H'j-h'\ =0(u(h';l/J)). 

Jt 

Let T be a large real number. Applying integration by parts to the left side of (j2.2j) and 
using (jl.4j) and (jl.5|) . we obtain 

hd^ T = f h{y) dM{y; T) 

T JT 



(3.7) 



Hj(y) dM(y; T) + / (%) - Hj(y)) dM(y; T) 
Jt 

Hj(y)dM(y;T)- [ M(y;T) (h'(y) - H'j(y)) dy 
Jt 

f E A{ P c * E *7 + (W) lo s T / \ h ' - H 'A > 



l<|i|<2J-2 0<7<T 



where we have written Xj = e 271 "- 7 ". 

To apply Lemma we write x 11 = x p ~ 1 ^ 2 + x n (l — x' 3-1 / 2 ). Let 5 = 50 lo f^°^ T ■ Assume 
that 

< 5losx < 1. 
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Note that if (3 + 27 is a zero, then 1 — (3 + ry is also a zero. We obtain from Lemma El the 
bounds 



and 



Thus 
(3.8) 



x^il-x^ 1 ' 2 )^ 



x 



(3-1/2 



0<7<T 

\P-h>s 



0<7<T 



<x 5 A^(l/2 + 5,T) + logx I 

Jin 

T- s/4 + \ogx J 



< TlogT 
T 

< 



x ff -V 2 N(a,T) da 

1/2+5 
1/2 /rpl/4" 



X 



(10 



+ 2-1/2-/3 



2) 



0<7<T 
|/3-||<5 



0<7<T 
|</3<i+5 



«log 2 x ^ (/3-1/2) 5 



0<7<T 
|</3<i+5 

-5 



<log 2 x / 9N(l/2 + 6,T)d6 
Jo 

T log 2 x 



logT 



{x h -x p ~ 1/2 ) < 

0<7<T 



T log 2 x 
logT 



We apply these estimates with x = Xj, 1 < j < 2 J — 2, noting that ar!2j = a; - 17 and c_j = Cj. 
By Lemma ^ we obtain 



hdn a T = 3? 

7T 



(3.9) 



2 ^ A { pc 3 K(n X] ) sin (Tlog^/r^)) 
^ v^i Tlog^/n^.) 

2J-2 2J-2 



+ ( ^ E i c ih 1/2 lo g 2 ( 2 ^ T ) + E m lo § 2 x 



i=i 



+ o(V(T) logTjf \h'-H'j\\. 
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Now put 

(3.10) J= UlogT^ + ij. 

Using (|3.5j) . (|3.10j) and \ogXj = O(j), we find that the first error term in (|3.9|) is 

J 2 



O 



X 2J-2 lo § 2 T 



T 



+ 



logT 



o(l) (T-oo). 



We now estimate the second error term in (|3.9|) . Since hi E Li(T), j T \h f — H'j 
J -> oo. Thus, if D*(r) < 1/logT, then 



as 



(3.11) 



£>*(T)logT / 1^-^1 = 0(1; 
./t 



Similarly, (EHTJ) follows if (JHU holds by (J3~K|) . 

Consider now the main term in (|3.9|) . This sum is absolutely and uniformly convergent in 
T and each term with n Xj 7^ Xj tends to zero as T — > cxo. Thus, the aggregate of such terms 
is o(l) as T — >• 00. Therefore, by ()3.8|) and (|3.11|) . 



f 1 

/ /l(iyU a t = 3? 



E 



+ o(l) (T^oo) 



In the case where a is not of the form for positive integers a, g and prime p, A(xj) = 
for all j and this finishes the proof. If a does have this form, A(xj) = unless q\j, in which 



case A(xj) = log p. Thus 



logp 



7T 



E 



/1 (J) r 



m<(2J-2)/q 



,am/2 



+ o(l) (T^oo) 



It is easily seen that (since h and g a are absolutely continuous) that the right side approaches 
/ hg a as T — > 00. This completes the proof of Theorem 1. 



4. Proof of Theorem [7] 

First, if CL\ , (32, ... IS not uniformly distributed in [0, 1), by Weyl's criterion the conclusion 
of Theorem 6 holds with h(y) = e 2mky for some integer k, regardless of the function /. 

From now on assume that CLiy CL2, ... IS uniformly distributed in [0, 1). Let r(x) be a 
C°°(1R) function such that r(x) = for \x\ > 1/2, r(0) = 1 and r is monotone on [—1/2,0] 
and [0, 1/2]. For any real number u, and any v > 0, 5 > 0, denote by r UjVy $ the function 
defined by 
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Thus r UjVj g is nonnegative, is C°°(IR), it is supported on an interval of length 8 centered at 
u, and one has r U)V; g(u) = v. Note also that 



/ 

7 k 



r u ,v,s(x)dx — v 8 / r(x)dx 



We now proceed to construct a function h as in the statement of the theorem. It is enough 
to define h on the interval (t — l,t]. We set h{t) = 0. Next, we write the interval (t — l,t) 
disjoint union of intervals 



1 1 



, k — 1, 2, 3, ... , 



and we define h inductively on each of the intervals Ik- We also set 

Jk 



7 5 



2*;+2' 2 fc+2 

for any k > 1. Thus is contained in and length (</&) = l/2 fe+1 . For each > 1, we will 
select a positive integer nk, and two real numbers Vk > 0, 5k > 0. Then /i will be defined on 

h by 

Here the sum is restricted to distinct values of a n . Thus if a n is the same for several values 
of n < rik, only one of these values of n is taken on the right side of (|4.1j) . We will choose 8k 
to be smaller than l/2 fc+2 . Then the right side of (|4.1j) will indeed be supported inside the 
interval Ik- Moreover, h will vanish on a small interval around each of the points t — ^ with 
k > 1. Note also that on each Ik, h is a finite sum of C°° functions, so after this construction 
is complete h will be C°° on T\ {£}. For each k, we only choose 8k after rik has already been 
chosen. Then we let 8k to be small enough so that \a n — a n >\ > 8k for any n,n' < rik with 
0"n' 7^ a n- This will make the supports of the functions r aniVki $ k from the right side of (|4.1|) to 
be disjoint. As a consequence, one will have < h(x) < Vk for any x G Ifc. If we choose the 
sequence {vk)k>i to be decreasing to 0, the function h will be continuous at t. If furthermore 
the sequence {vk)k>i is chosen so that 

lim 2 k v k = 0, 

then h will be differentiable at t, and h'(t) = 0. We put Vk = l/3 fe , so the above condition 
holds. It remains to construct the sequences (rik)k>i, and (8k)k>i- Fix a k and assume 
that rij and 5j have already been defined for j = 1, . . . ,k — 1. This means that ft, has been 
constructed on the interval Ej. := I\ U • • • U Ifc-i- Let n' k be such that for any n > n' k we 
have f(n) < l/7 k . Now, since the sequence (a n ) n >i is uniformly distributed, there exists an 
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n'l such that for any n > n'l one has 



Kj<n Jtj k 



i<i 



< 



*jk ' 



The fact that (a n ) n >i is uniformly distributed also implies the existence of an n!' k ' such that 
for any n > n'l', 



#{m <n : a m e J k } > 



n 



2 k+2 ■ 



We put rik = max{n' fc , n" k , n'l!}. Then by the above it follows that 

f(n k ) < , 



i<i<« fe 



1 

— 'jk 



and 



#{j < n k : Oj G J fc } 



i<i<«fe 

CLj£l k 



n k 



Vk > 



2k+2^k ' 



Lastly, we choose 5k > to be small enough so that it satisfies all the requirements stated 
so far, and such that we also have 

1 



h 7 

This completes the construction of the three sequences (njfe)fc>i, (v k ) k >i, (5 k ) k >i, and thus 
also the construction of h. It remains to check that the inequality from the statement of 
the theorem holds for infinitely many n. Take n = nk, and break the interval (t — l,t] as a 
disjoint union of three intervals, (t - 1, t] = E k U I k U T k , where T k — (t — l/2 k , t\. We also 
break accordingly the sum and the integral from (|2.3|) . 



1 n k 

-E>( 



and respectively 



We know that 



h+ I h+ I h. 

Ek J Ik 



Z Ek - J Ek h < l/7 k , E 4 > 1/(4 • Q k ), < J h h < l/7 k and f(n k ) < l/7 k . 
Also S^ fc > 0, and since for any m > k we have < f T h < 1 /7 m , it follows that J T h < 2/7 k . 
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We conclude that for large k, 

^ !>,)-//> >/(„,,, 

which completes the proof of the theorem. 
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